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Abstract 

. The general SU (N) form factor formula is constructed. Exact form factors for 

QJ I the field, the energy momentum and the current operators are derived and compared 

with the 1/A^-expansion of the chiral Gross-Neveu model and full agreement is found. 
As an application of the form factor approach the equal time commutation rules of 
arbitrary local fields are derived and in general anyonic behavior is found. 
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^ '. 1 Introduction 
o 

Quantum chromo dynamics, the theory of the strong interactions, is a non abelian gauge 
theory based on the gauge group SU{3). It was first pointed out by 't Hooft [U [2] that 
many features of QCD can be understood by studying a gauge theory based on the gauge 
5^ I group SU (N) in the limit N oo. One might think that letting N oo would make the 
analysis more complicated because of the larger gauge group and consequently increase 
in the number of dynamical degrees of freedom. Also one can think that SU{N) gauge 
theory has very little to do with QCD because ^ oo is not close to = 3. However it 
is well known that the expansion provides good results which can be compared with 
experiments [5]. 
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One of the most important trends in theoretical physics in the last decades is the 
development of exact methods which are completely different from perturbation theory. 
Resolution of the strong coupling problem would give us a full understanding of the struc- 
ture of interactions in nonabelian gauge theory. One promising possibility of overcoming 
the limitations of perturbation theory is the application of exact integrability. From this 
point of view the two dimensional integrable quantum field theories are in a sense a lab- 
oratory for investigations of those properties of quantum field theories, which cannot be 
described via perturbation theory. 

The chiral SU (N) Gross- Neveu [3] model given by the Lagrangian 

is an interesting 1-1-1 dimensional field theory that can be studied using the 1/iV expansion. 
The model is asymptotically free with a spontaneously broken chiral symmetry, and so 
shares some dynamical features with QCD. Gross and Neveu [1] investigated the model 
using an 1/A^ expansion. Apparently a chiral ?7(l)-symmetry is spontaneously broken, 
the fermions acquire mass and a Goldstone boson seems to appear. This is of course 
not possible in two space-time dimensions and severe infrared divergences appear due the 
"would-be-Goldstone boson" . However, it has been argued by Witten [5] that dynamical 
mass generation can be reconciled with the absence of spontaneous symmetry breaking. 
There exist further (different) approaches to overcome these problems and to formulate a 
1/N expansion P, [3, |B] (see also [9J). On shell they all agree and are consistent with the 
exact S-matrix ([2]) . We follow here the approach of Swieca et al. [S] where additional fields 
are introduced in order to compensate the infrared divergences. The authors claim that 
since the physical fermions have lost not only the chiral U{1) but also the charge U{1) 
symmetry, they transform accordingly to pure SU{N). They propose an interpretation 
of the antiparticles as a bound state of iV — 1 particles. Furthermore this means that the 
particles satisfy neither Fermi nor Bose statistics, but rather carry "spin" s = | (1 — 1/A^). 
As a consequence there are unusual crossing relations and Klein factors. 

In this article we will focus on the SU{N) Gross-Neveu form factors using the "boot- 
strap program" fiUi [TT] . We provide here some examples, calculate the form factors 
exactly and compare the results with field theoretical 1/N expansions. We emphasize 
that in addition to the operators in the vacuum sector, such as the energy momentum 
tensor and the current, we also consider anyonic operators as the fundamental fields. We 
also derive the equal time commutation rules for local operators which are in particular 
complicated due to the unusual crossing formulae related to the Klein factors. 

The general form factor of an operator 0{x) for n-particles, which is a co- vector valued 
function and can be written as [12] 

fs{o)=k^{o) n ^(%) 

l<i<j<n 

where 9 = {9i, . . . ,9n) is the set of rapidities of the particles a = (ai, . . . , a„). The scalar 
function F{9) is the minimal form factor function and the K-function K^[9_) contains 
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the entire pole structure and its symmetry is determined by the form factor equations 
(i) to (v) [13]. To construct the K- function we must apply the nested off-shell Bethe 
ansatz to capture the vectorial content of the form factors. This solves the missing link 
of Smirnov's [14j formula for the SU {N) form factors, where the vectors were given by an 
"indirect definition" characterized by necessary properties but not provided explicitly. We 
note that SU (N) form factors were also calculated in [TH US, HB] using other techniques, 
see also the related papers [TTl [IB] . Our results apply not only to chargeless operators 
such as the energy momentum and the current operators but also to more general ones 
with anyonic behavior. We believe that our integral representation, besides of being 
appropriate for a comparison with field theoretical expansions, may also shed some 
light for a better understanding on the correlation functions of models with more general 
(anyonic) statistics. 

The paper is organized as follows: In section [2] we present the general setting concern- 
ing the SU{N) S-matrix, the nested off-shell Bethe ansatz and the chiral Gross-Neveu 
Lagrangian field theory. We review known results and derive some further formulae which 
we need in the following. In section [3] we construct the general form factor formula and 
present some examples in detail, such as the form factors of the energy-momentum tensor, 
the Dirac field and the SU{N) current. In section H] we compare our exact results against 
perturbation theory of the SU{N) Gross-Neveu model. In section [5] we present the 
commutation rules of the fields. Our conclusions are stated in section [61 In appendix [Xj 
we provide the general proof of the bound state form factor formula and in appendix [B] 
the commutation rule of two fields (in general anyonic) is proved. 

2 General setting 

The particle spectrum of the chiral SU{N) Gross-Neveu model consists of A^ — 1 multiplets 
of particles of mass rrir = mi sin {rn/N) / sm{n/N), which correspond to all fundamen- 
tal SU{N) representations of rank r = 1, . . . , — 1 with representation spaces V^*^^-* of 
dimension (^). Let (a) = (ai, . . . , a^), (1 < ai < ■ • ■ < ctr < A^) be a particle of rank r. 
We write 

N-l 

[a) eV = ^V'-'\ ~c(^) 

r=l 

where the (a) form a basis of V. A particle of rank r is a bound state of r particles 
of rank 1. The antiparticle corresponding to (a) is (a) = {ai, . . . ,aN-r), (1 < «i < 
■ ■ ■ < aN-r < N) (of rank N — r) such that the union of the set of indices satisfies 
{«!, . . . , a J U {«!, . . . , aN-r} = {1, . . . , N}. 

2.1 The S-matrix 

The S-matrix for the scattering of two particles a, P (of rank 1) [191 [201 [HI [8] is 

si}{9) = 5i5',b{e) + 5i5;cie) (2) 
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where 6 = 6i — 62 is the rapidity difference and = fn (cosh6'i 2, sinh^i 2). The amph- 
tudes satisfy 

r fi - ^) r Ci - ^ + ^) 
aie) = bie) + 0(6) = ■ 

We also need the S-matrix for the scattering of a bound state (p) = (pi, . . . , Pn~i), (1 < 
Pi < ■ ■ ■ < ^ N) (of rank — 1) and a particle a (of rank 1) 

^SlV) = (-1)"^-' (^S^^ Km -6) + C^('^)C(,). c(vr. - 6)) (3) 
where the charge conjugation matrices are defined by 



Qai(a2...ojv) _ (~((ai...ajv_i)«jv — ^_]^^Af-lgaia2 



JV 

...Oiv 



with eaj...Q,^ and e"^' '"^ totally anti-symmetric and ei...Ar = e^ --^ = 1. Formula ([3]) is 
obtained by applying iteratively the bound state fusion method [21j to d^D- 

For later convenience we extract the factors a{6) and {—l)^~^b{i7r — 6), respectively 

st}{e) = a{9)S%{e) (4) 
W = {-lr'%^'r^ - 0)&e) (5) 



such that 



S%{e) = 6l6lm + 6i6}c{e) (6) 
Si!>)lH = + C^('^)C(,)«d(^) (7) 

6 — 17] d — 17] 0{77T — UJj m — UJ N 

where S^^^ = 6p^6p^...6p^zl ■ Below we will also use for the matrices (jH]) and ([7]) the 
notations 5'i2(6') and Si2{9), respectively. 



2.2 Nested "ofF-shell" Bethe ansatz 

The "off-shell" Bethe ansatz is used to construct vector valued functions which have 
symmetry properties according to a representation of the permutation group generated 
by a factorizing S-matrix. In addition they satisfy matrix differential [22] or difference 
[23] equations. For the application to form factors we use the co- vector version ^ {6_) G 

K...n = (OLi (^. ec, ^ = i,...,n) 

K,,,ij„X- ■ ■ ■ ■ ■) = -f^...ji...(- • • , Oj, Oi, ■ ■ ■) Sij{9ij) 
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where 9^ = {9i, . . . , 9i + 27ii, . . . , 6'„) (see below and e.g. [231 [IS])- We write the components 
of the co-vector K-^ ^ as where a = ((an, . . . , airj, • • • , (ctni, • • • , «nT-„)) is a state of 
n particles of rank ri, . . . , r„. 

The nested SU{N) "off-shell" Bethe ansatz for particles of rank 1 has been constructed 
in [13]. Here we need a more general case. 



Nested "off-shell" Bethe ansatz for particles of rank 1 and — 1 : We consider 
a state with n particles of rank 1 and n particles of rank — 1 and write the off-shell 
Bethe ansatz co-vector valued function as 




where a = (ai, . 
9 = {9i, . . . , 9n) 



•,an), (p) 



((pi), . . . , (pn)) = ((Pll, . . . , PlJV-l), . . . , {Pnl, PnN-l)), 



and z 



This ansatz transforms the 



complicated matrix equations to simple equations for the scalar function k{9,ui,z) (see 
[13] and below). The integration contour Cqi^ can in general be characterized as follows: 
there is a finite number of complex numbers ai{9_), bj{9_) such that the positions of all 
poles of the integrand are of the form 



(1) : ai{9) + 2Tiik , keN 

(2) : bj li) - 2TTil , leN 



(9) 



and Cqui runs from — oo to +oo such that all poles (1) are above and all poles (2) are below 
the contour. This contour is just the same as the one used for the definition of Meijer's 
G-function. It will turn out that for the examples considered below the form factors can 
be expressed in terms of Meijer's G-functions. 

The state ^E'a(p) in ([8]) is a linear combination of the basic Bethe ansatz co-vectors 



'^aip){9,ui,z) 



Lpia)_{z,ui)^^( 



,ui,z), with 1 < Pi, aij = 1 . 



(10) 



As usual in the context of the algebraic Bethe ansatz 
co-vectors are obtained from the monodromy matrix 



.71,1. ..71, 



fl{9,Ul, 9o) = Sio{9i — 9o) ■ ■ ■ Snoi9n — ^^o)'S'io(co'i — ^^o) 
^l...n,l...n(^) ^) z) Bi n l n f^{9, U, z) 



C: T .(t/,u,z 

1.. .71,1.. .71 ^— ' ' 



D, J _ At^,uj,z} 

l...n,l...n,f)^—' — ■ 



the basic Bethe ansatz 

■ Sno{u!n — 9q) 

, 2<(3,(3'<N. 



where the S-matrices Siq and Sjo are given by ([6]) and ([7j). As usual the Yang-Baxter 
algebra relation for the S-matrix yields the typical TTS-relation which implies the basic 
algebraic properties of the sub-matrices A, B, C, D. 

Here not only one reference co- vector exists. The space of reference co- vectors, defined 
as usual by 

n W 5^ = , 
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is (A^ — 1)" dimensional and is spanned by the co- vectors for all (a) = ((cth, . . . , ai^^i), 

. . . ,{ani, • • • , CTfiN-i)) with cTji = 1 < crj2 < ■ ■ ■ < cTjAf-i < N. They are eigenstates of A 
and Dp 

Q^^A{e,u,z) = 

i=l 

where the indices 1 . . . ra, 1 . . . n are suppressed. The basic Bethe ansatz co- vectors in ffTOl) 
are defined as 

^, ^) = h^^c^- a, u, zj--- {e^u^z,)) (11) 

where 1 < /3i < A^. 

The technique of the 'nested Bethe ansatz' means that for the coefficients Lj3(^„-j{z,ui) 
in (II up one makes the analogous construction as for Ka{p){6.,uj_) where now the indices 
P, (cr) take only the values 2 < Pi < N and an = 1 < < ■ ■ ■ < CTi^^i < N. This 
nesting is repeated until the space of the coefficients becomes one dimensional. It is well 
known (see [23]) that the 'off-shell' Bethe ansatz states are highest weight states if they 
satisfy a certain matrix difference equation. If there are only n particles of rank 1, then 
the SU{N) weights are 

w = {n - ni,ni - n2, . . . , 7x^-2 - nN-i,nN-i) (12) 

where rii = m,n2, . . . are the numbers of C operators in the various levels of the nesting. 
If in addition there are n particles of rank — 1 the SU{N) weights are 

w = (n - ni, ni - ^2, . . . , nN_2 - Un-i, n^^i - n) + n{l, . . . , 1) (13) 

because — 1 particles of rank 1 yield a bound state of rank A^ — 1 and at the l*^ level 
the number of C operators is reduced by A^ — / — 1 (see appendix [XI). 



2.3 Minimal form factors and 0- function 

To construct the form factors we need the form factor functions F {6), G {$) and the 
function 0(^). The form factor functions F (6) and G (6) for two particles of rank 1 and 
for one particle of rank 1 and one of rank A^ — 1, respectively are 





oo 
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They are the minimal solutions of the equations 

F{9) = F{-e)a{e), F{m -9) = F{in + 9) 
G{9) = -G{-9)b{7it - 9) , ^(m -9) = G^m + 9) 

where a{9) and hijii — 9) are the highest weight amplitudes of the corresponding channels 
of ([2]) and da]). The 0-function 

is a solution of 

N-2 N-1 

JJ (-0 - kir]) Y[F{9 + kiri) = 1 

k=0 k=0 

which follows from the assumption that the antiparticle of a fundamental particle is a 
bound state of — 1 of them (see below and [IS])- The constants c and c' in (fn|) and 
follow from ([TBI) and (fTTl). 



(17) 



2.4 Chiral Gross-Neveu model 

Swieca et al. p] wrote the fermion fields ipi{x) in the Lagrangian ([T]) in bosonic form. In 
order to extract the real particle content of the theory, they introduced in addition the 
"physical" fields 



'4^j{x)ICi for 2 = j 
-ijj{x)ICi for 2 7^ j . 



where V5j(x) are free canonical zero-mass fields and /C, are Klein factors satisfying 

JCiiJjj (x) 

The fields ip satisfy (with a suitable normal product prescription A/") 

1 . . ^ 

such that 



K,i}j{x) = {-\Y-^i)j{x)K,. (19) 
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Equation (IT5]) means that antiparticles should be identified with bound state of — 1 
particles and the creation operators of the antiparticle 6]^ and of the bound state a|^^ are 
related by 

lCailC = {-lf-'ai. 
The "physical" fields satisfy the anyonic commutation relations 

Ux)Uy) = Uy)Uxy^^'''^"'~'''\ for - y)' < 

with "spin" s = | (1 — l/A^)- This implies that the "physical" S-matrix is related to the 
one of ([2]) by P [26] 

As a consequence the abnormal crossing relation ([3]) transforms to a normal one. The 
bound state S-matrix satisfies 

Sl^liO) = {-lf-'C^,),St{'r^^ - 0)C'^'^''\ (20) 
Therefore the physical crossing relation is 

with Caa' = ^aa'i C'^'^ = ^'^''^ ■ 

3 Form factors 

For a state of n particles of rank ri, . . . , r„ with rapidities 9_ = (9i, . . . , On) and a local 
operator 0{x) we define the associated form factor functions F^{6_) by 

( I 0{x) I ^1, ... , ft„ = e— (^^1+-+^")^^ (^) , for > ■ • ■ > 9^. 

where again a = ((an, . . . , air J, ■ ■ ■ , (a„i, . . . , a-nrn))- For other arrangements of the 
rapidities the functions F^{6_) are given by analytic continuation. The co- vector valued 
function F^{6) satisfies the form factor equations (i) - (v) (see [121 [HI [271 [2HI [13]) and 
can be written as [12] 

F°{e) = K^{e) n (21) 

l<i<j<n 

where Fr^r {&) are the minimal form factor functions. For particles of rank 1 and — 1 
they are given by Fn(^) = F^.i^-iiO) = F{e) and F^_ii(^) = Fi^-i(^) = G'(^) of m 
and (1151) . respectively. In [13] the form factors of the fundamental particles of rank 1 have 
been constructed. We shortly recall the results. All other form factors can be obtained 
from these by applying the bound state fusion procedure which is given by the form factor 
equation (iv) (see e.g. [13]). 



3 FORM FACTORS 



9 



Form factors for particles of rank 1: The K-function in f l21l) is given by the nested 
"off-shell" Bethe ansatz ([HD for the special case n = and a special choice of the scalar 
function k{6_,z_) such that the form factor equations (i) - (v) are satisfied 



with 



K. 



o, 



^ I dzi- 



Ce 



dZmh{e,z)p''{e,z) -^aiO^z) 



Ce 



(22) 



i=i j=i 

r{z) = ^ 



n 

l<j<j<m 



T[Zi - Zj^ 



(23) 



z)<p{-z) 



The integration contour Ce is defined by ([9]). The dependence on the operator O enters 
only through the p-function p'~'{6_,z) which has to satisfy simple equations (see [291 EDI 
EH [IS])- The K-function is in general a linear combination of the fundamental building 
blocks [2Sl[SniEI] given by ([22D- Here we consider only these cases where the sum consists 
only of one term. 



The p-function: The co- vector valued function \l/a(^, z) is expressed as in (ITOj) for 
n = by an Lp{z) for which the nesting procedure is applied. The final form is (up to a 
constant) 

Faii) = n^(%) / ^^^'^ • • • / d^^''''^KO.ip''{0-.i) ^M^i) (24) 

where z = z^^\ . . . ,z^^~^^- The Bethe state ^a{^,z) is obtained by the nesting procedure 
(see (Uni), (HH) and p]) 

^a(^,i) = ^''-'\z^''-'\z^''-'^) . . . ^f'Hz^'\z^'^Mie,Z^'^) . 

In general the p-function (see [13]) depends on the rapidities 9 and all integration variables 
z^''\ Let the operator 0{x) transform as a highest weight SU{N) representation with 
highest weight vector 

Because of SU{N) invariance the weight vector of the co- vector F^{6_) is then 

w={wf,...,w'^)+L{l,...,l) (25) 
= (ra - rii, ni - 722, . . . , - tiat-i, tiat-i) 

where f lT2|) and the fact, that the weight vector (!,...,!) correspond to the vacuum sector, 
has been used. In [13] was shown that the p-function has to satisfy a set of equations 
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in order that the form factor (1241) satisfies the form factor equations. In particular to 
guarantee the transformation properties of the operator the following periodicity relations 
have to be valid 

= (-l)-'?+<ipO(e,...,zf^ + 2m,...) ^ ^ 

where 

= (-l)(^^i)&"'?/^]-^-2-f 
af = e'-(i-i/^)«'' . (27) 

The charge of the operator O is defined by Q'-^ = n mod N and af is the statistics factor 
of O with respect to the fundamental particle of rank 1. The sign factors /af = ±1 
and (— 1)""' +"^'+1 = ±1 in (12^ follow [13] from the sign (— 1)(^~^) in the unusual crossing 
relation fl20|) (related to the Klein factors of f|T9l) ). 



3.1 General form factors of particles of rank 1 and — 1: 



Using the bound state procedure (see appendix which means taking residues of fl^T]) or 
( 122|) one derives the form factors and K-functions for n particles of rank 1 with rapidities 
6_ and n particles of rank — 1 with rapidities u. As usual we split off the minimal part 



l<i<j<n i=l j=l 

The K-function is given by a nested 'off-shell' Bethe ansatz 



n 

l<j<j<n 



{2i 



dzi ■ 



dzm h{9, z) {9, u, z) ^fo(p) {9, lj, z) 



(29) 



where h{9_,z) is the scalar function fl2^ . Note that this h-function does not depend on a;. 
For the SU{N) S-matrix the function (j){9) is given by (|T6|) . The integration contour C^^^ 
(see Fig. [1]) has been defined in the context of (jH]). 



Nesting: The state ^I/a(p) in (|29|) is a linear combination of the basic Bethe ansatz 
co-vectors ffTTj) 

^a{pii9,ui,z) = Lp^{z,ui)^^{9,ui,z) , with 1 < pi, aij = 1 

where -C//3(cr)(2;, cj) satisfies again a representation like (l29l) . This nesting is iterated until 
all Pi = N and all {a)i = {1,2, . . . , N — 1). Only for the highest level Bethe ansatz the 
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uj + Sir 



'6*2 + 27ril 



.0i + 2vrz(l-i) 



u + in — 



uj — i-K — 2nij^ 



02 - 27rzi 



* 1 - 27ra 



00 - 2ni 



Figure 1: The integration contour Cg-^g^^ for two particles and one bound state. 



h-function depends on u. The final result is 

N-2 n rijv-i 

1=0 1=1 j=l 

= KMa) " G - i - 1?) ■ 

The complete Bethe ansatz state is 

where (r^) denotes n highest weight bound states (r^ji, • • • , fliN-i) = (1, 2, . . . , — 1). The 
p-functions in ( 129|) and (130|) are obtained again by the bound state procedure from a 
solution of (126|) for n = 0. In particular for n = 1 (with the replacements in (l26l) 0_-^ 6_,(f 
and z(') where ^ = (v?i, • • -^ipN-i) , y^'^ = (vi, ■ ■ .,yN~i~i) , / = 1, . . . , - 2) 

Here ?/f ■* = \ i = 1, . . . , N — 1 — I and (/^^ = + fcir^ — iir. The proofs of the statements 
of this subsection and more details can be found in appendix \M 

3.2 Examples 

To illustrate our general results we present some simple examples. In addition, we also 
derive the 1/N expansion of exact form factors for the purpose of later comparison with 
the 1/iV-perturbation theory of the chiral SU{N) Gross- Neveu model. 
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The energy momentum tensor: For the local operator 0{x) = Tp^{x) (where p,cr = 
± denote the light cone components) the p-function for n particles of rank 1 (as for the 
sine-Gordon model in [28]) 

n m 
i=l i=l 

satisfies the equations (126|) with w'^ = (0, 0, . . . , 0). For the n = N particle form factor 
the weight vector isw = (l,l,...,l,l). Due to (fT2l) there are ni = N — I integrations in 
the l-th level of the off-shell Bethe ansatz. 

We calculate the form factor of the particle a and the bound state (A) = (Ai, . . . , Aat-i) 
of A^ — 1 particles. We apply the bound state formulae fl28l) and fl29l) for n = n = 1. Due to 
f fT3|) there is just one integration in every level of the nested Bethe ansatz [1 = 1,..., N—1) 



FH:^{e,^) = Kl[:^{e,u)G{e-u;) 

K^:^{e,u;)=Nr (e^' + e^-) [ dz4>{9 - z)e''^^^(^,){e,u;, z) (32) 

JCg_ 

(9, CO, z) = L« ) {z, (9, u, z) 



where G{9) defined in f[T^ is the minimal form factor function of two particles of rank 1 
and N — 1. The integration in every level of the nested Bethe ansatz (/ = A^ — 2,...,1) 
can be solved iteratively 

= 2^) with/?>/, (/i) = (1,2,...,/,*,...,*) 

..V-,^c,r(i.^)r(-i.l-^). (33) 

The remaining integral in (1521) may be performed (see appendix [X] ) with the resu 110 

( I T^'^(O) I 9, u Y:^,^ = 4m^6.(A)ei(^+-)(^+-+-) '^^^j ~ T G{9 - u) . (34) 
^ ' u — uj — m 

Similar as in [28] one can prove the eigenvalue equation 



j dxT^\x) - Y,pf^ I ^^1, . . . , = 



for arbitrary states. 



'^In [3H[53] this result has been obtained using Jackson type integrals. 
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The fields ipaix)'- Because the Bethe ansatz yields highest weight states we obtain the 
matrix elements of the spinor field ip{x) = ipi{x). The p-function for the local operator 
ijj^'^^x) for n particles of rank 1 (see also j27j ) 

/^'(^,.)=exp±i (f:^.-(i-^)j:e. 

\i=l ^ ^ 1=1 

satisfies the equations (l26l) with = (1, 0, . . . , 0). For example the 1-particle form factor 
is 

The last formula is consistent with the proposal of Swieca et al. [HI E] that the statistics 
of the fundamental particles in the chiral SU{N) Gross- Neveu model should be cr = 
exp (27ris), where s = | (l — is the spin (see also (127|) ). For the n = N + 1 particle 
form factor there are again ni = N — I integrations in the /-th level of the off-shell Bethe 
ansatz and the SU{N) weights are w = (2, 1, . . . , 1, 1). Due to (IT^ there is again just 
one integration in every level of the nested Bethe ansatz. Similar as above one obtains 
the two-particle and one-bound state form factor 

Fi'Uo) = Kt;^4e_)F{e,,)G{er,)G{e,,) 

KfftiUo) = iV^e^K(i-^)(^^+^^)+^^3) f ^^^(^^ _ ^)^(^^ _ ^)e±^^^,^(,)(£ z) (35) 

^a/3(A)(^, Z) = L^if.){z, e3)^l^p%^{e, Z) , with 1< 7, Ai = 1 

where the function L^(^^{z^9^) = e^(^f^)L^^\6^ — z) is the same as in ([3SD above. We were 
not able to perform this integration, however, the result can be expressed in terms of 
Meijer's G-functions 



^ _L 1 ^ _L 1 ^ I 3 _ J_ 
27r« Af' 27ri TV ' 27ri 2 



and K2'' ' is obtained by the form factor equation (i). 

1 /N expansion of the exact form factor: We consider the connected part of the 
matrix element 

^ e, I ^t\x) I e, )™"- = c(^)-F^;(^)(^i, e,, e, - r-n) . 

Instead of the field ip we consider the operator Os = {—i {i'-fd — m) ip)^ 
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For N oo we expand the minimal form factors 

F{e) = — sinh + 0{1/N) , G{e) = l + 0{1/N) , 

71 

perform the integration in (l35i) and obtain (after a lengthy calculation) 

^0 _ 2mi7r sinh 6^13 / 1 5 ^ \ fa \ niAT-'iy 



We use the following conventions for the 7-matrices 



-1 

1 



and for the spinors 



m(p) = V^( ) , V{p) = ^l\ _^g/2 )• (37) 



In section m below we compare this result with the 1/A^-expansion of the chiral SU{N) 
Gross- Neveu model in terms of Feynman graphs. 

The current J^pi^)'- The SU{N) current J^i^p-^ix) transforms as the adjoint representa- 
tion with the weig ht vector = (2, 1, . . . , 1, 0). A gain, because the Bethe ansatz yields 
highest weight states we obtain the matrix elements of the highest weight component 

where we have introduced the pseudo-potential ip{x). We start from 

with z = z^^\ . . . ,z^^~^^- The proposal for the p-function for n particles of rank 1 (see 
also [27J) 



ni "-JV-1 



( Y: exp «. exp 1 5: 9, - 5: ' - .f'' 

,1=1 / \j=l i=l i=l 



satisfies the equations (l26i) with w'^ = . 

We calculate the form factor of the particle a and the bound state (A) = (Ai, . . . , Aat-i) 
of — 1 particles with weight vector w = (2, 1, . . . , 1, 0). We apply the bound state 
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formulae (1281) and fl29l) for n = n = 1. Due to fll3p there is no integration in each level of 
the nested Bethe ansatz (/ = 1, . . . , — 1) and 

The form factor for the SU (N) current is therefore 

Kll^ (0, ^) = ( I J,t)(0) I 9, u = ±N,6i5^^ (e±^ + e^^) ^^G{e - u) 

= 5f 5;^; v{u;hM0) G{e - u)/G{zrr) . (38) 

Also here we calculate the 1 / iV-expansion of the exact form factor for later comparison 
with the 1/A^-perturbation theory of the chiral SU{N) Gross-Neveu model. Using the 
expansion of the minimal form factor function 

obtain the 1/A^ expansion of the exact the SU{N) current form factor as 
fJ;' [6, u) = {0 I J|(^)(0) I e, u = 5f 5g v{u:)^^u{e) G{e - u)/G{zn 

4 The chiral SU{N) Gross-Neveu model 



we 



0{N 



~2\ 



Let the fermi fields ipaix), (a = 1, . . . , A^) form an S't/(A^)-multiplet. The field theory is 
defined by the Lagrangian [1] 

or equivalently 

C{ip, tjj, a, tt) = ^{i'yd — a — i'j^Tc)^/; — ^g^^(a^ + tt^) 
where o"(x) is scalar and 7r(x) a pseudoscalar field. The field equations for these fields are 



4 THE CHIRAL SU{N) GROSS-NEVEU MODEL 



16 



4.1 The 1/N perturbation theory 

Using the bootstrap program and the results of [19], the S-matrix i.e. the on-shell solution 
of the model has been proposed in [3, [8]. It is well known [H El [3, [8] that the naive 1/A^- 
expansion of the chiral Gross-Neveu model suffers on severe infrared problems. In [7, 8J 
two different approaches to overcome these problems were proposed and it was shown 
that the exact S-matrix was consistent with both. We will show that an off-shell quantity 
as our solution for the three particle form factor of the field ip{x) is also consistent with 
the 1/A^-expansion of |8J- Without presenting details we note that we do not obtain 
consistency with the approach of [7]. Since in the hterature (see e.g. [8l[33l[9j) there are 
some errors and misprints we present a detailed derivation of the approach of Swieca et 
al. 

The generation functional of Greens's functions for the chiral Gross-Neveu model is 

Z(a) ^ / # * exp , (AU, ^, .) + » + «) (39) 

with the action Ai^ip, i^, cr? ^r) = J d?x C^ip, V", o", ^r). In eq. (139|) and in the following we use 
a short notation of the x-integrations e.g. ^ip = J (fx^{x)xl'{x). 

When quantizing the model, severe infrared divergences appear due to the "would-be 
Goldstone boson" vr. Following Kurak, Koberle and Swieca [8] we introduce two additional 
bosonic fields A{x) and B{x) quantized with negative norm. The A- field compensates 
the infrared divergences. In fact as we will see below that together with the infrared 
divergences of vr it decouples from the rest of the model. We replace the fermi fields by 

ipi^x) ip'{x) = exp i (^j^A^x) + B{x)) ip{x). 

The 5-field is introduced, in order not to change the statistics of the ^/'-fields. Finally we 
have the Lagrangian 

£ = tp'{i^d - /i)^' - y-^d'' + vr^) + lN{a-^ADA + p-^BDB) 
with n = a + i'j^Ti — 'j^jdA + 'jdB . 

The couplings a and (3 are unrenormalized, their renormalized values are a/tt. Performing 
the '?/''-integrations in the generation functional we obtain 

Z(^, f ) = J da d-K dA dB exp {iAeff{a, vr. A, B) - ^S^) 
with the fermi propagator S = i{i'yd — fi)~^ and the effective action 



Aeff{a,iT,A,B) 

= -iN Tr \n{i-fd - /i) - | [ d^x {g'^ {a^ + n^) - N{a-^ADA + (3~^BUB)) 
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The symbol Tr means the trace with respect to x-space and spinor space. The trace with 
respect to S'f/(A^)-isospin has been taken and gives the factor A^. We define the vertex 
functions F by 

-^effi^) = — / (fxi...(fXnT^^\xi,...,Xn)'fl{Xi)...LPniXn), 

n=0 ^' ^ ~ 

where (pi G {a — ctq, vr — ttq, A — Aq, B — Bq}. The values (Tq etc. are defined by the condition 
that Aeffif) is stationary at this point. This means that the one-point vertex functions 



x] 



6Aeff/Sa = etc. vanish 



r«(x) = N tT S{x,x) - g-'^ao = 

(x) = N ti {i-f^S{x, x)) - g-\o = 
t''}\x) = N ti (-7^7(95(x,a;)) + Na-'^DAo = 
r^J) (x) = N ti {-fdS{x, x)) + N^-'^UBq = 0. 

The three last equations mean ttq = Aq = Bq = and the first one implies = m with 

f (Pn i r, HL^ 

N / -— ^ tr g'^ao = ^ ctq = m = Me ^ , 



(27r)^ 7P — m 

where M is an UV-cutoff. There is the effect of mass generation and dimensional trans- 
mutation: the dimensionless coupling g is replaced by the mass m. The 1/iV-expansion 
is obtained by expanding the effective action at this stationary point. The resulting 
Feynman rules are given by the simple vertices 

V^{k) = i-t) , = 7' , VA{k) = 7V , Vnik) = -^k (40) 

and the propagators in momentum space 

X /,x "^^ 1 sinh(/) . 
A..(fc) = -- ,21, — , (41) 



N cosh^ i( 
in 1 



N sinh^ i0 V 



AAA{k) -- 

KA{k) = AAnik) = -2m- 



/ sinhi 



Nk"^ 

where k"^ = —Arn^ sinh^ To obtain the propagators one calculates the two point 

.(2) 



vertex functions F-- from the bubble graph of Fig 12] with the various vertices and uses 
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p + k 

p 

Figure 2: The bubble graph. 

A = zr(2)"\ In [8J it was argued that the unrenormahzed values of a and (3 are to be 
replaced hj a ^ oo and P ^/ti. In that limit the propagators are those of fHTj) . One 
observes that the A- and 5-propagators remain free and the infrared singularity in the 
7r-propagator disappears. 

As an example we consider the four point vertex function 

r^')?S(-P3, -P.,PUP2) = Si6jT^§{p,-ps) - Sl5f,T^^E{P3-Pi) 

where A,B,C,D are spinor indices, 7,5 are isospin indices and T is given by the 
Feynman graph of Fig. [31 Taking into account the contributions from all vertices (BOl) 





Pi P2 



Figure 3: The four point vertex 
and all the propagators fHTl) we obtain 

T{k) = Y,V^ik)A,,{k)V,{-k) (42) 

= -1 ® 1 A^^{k) + 7^ (8) 7^ KAk) - l^lk ® 7^7fc A^A(fc) 
- 7^ O 7^7^: A^a(A;) + 7^^ ® 7^ Aa^(A;) - 7A; ® jk ABB{k) . 
Inserting the expressions for the propagators we finally obtain 



in ( 1 sinh , , 1 / sinh (p 

+ ^ (7^7^; (g) 7^7^ + 2m7^ (g) 7^7 A; - 2m'j^'jk (g) 7^ + 7A; (g) 7/c) | (43) 

where the tensor product structure of the spinor matrices is obvious from Fig. [HJ We 
now apply these results to the examples of section [3] and investigate the three particle 
form factor of the fundamental fermi field and the two particle form factor of the SU{N) 
current in 1/A^-expansion in lowest nontrivial order. 
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The three particle form factor of the fundamental fermi field: For convenience 
we multiply the field with the Dirac operator, take 

Oosix) = {-i (i'jd - m) i'ix))^^ 

and define 

oJ(P3 I OosiO) I Pl,P2 = ^"^"^1/3(^12, ^^13, ^23) • 

By means of LSZ-techniques one can express the connected part in terms of the 4-point 
vertex function 

^connI/3(^12,^13,^23) = UciPs) T^Ba^i^P-i^ P3 ^ Pi ' P2,Pl,P2) Ua{Pi)ub{P2) ■ 

The lowest order contributions are given by the Feynman graphs of Fig. H] 

3 3 ^/-i 3 

/ 



0.(0)) = >--< - /\ + 
1 2 

Figure 4: The connected part of the three particle form factor of the fundamental fermi 
field in 1/N -expansion. 



FZnnlp = UciPs) {SasSfS^ T^sfe - Ps) " ^aySfss T^eiPs " Pi ) } (Pl fe) • 

where F is given by Fig. [3] and eq. and the spinor u{p) by eq. (^71) . It turns out that 
for pi, p2 and on-shell several terms vanish or cancel and we obtain up to order 

r) 7 2min f sinh 6*23/1 ^ 1 \ / \ 

which agrees with the result for the exact form factor fl5B]) . In [8J was shown that if the 
momentum Pa = pi + P2 — Ps is also on-shell then the expression (144|) is consistent with 
the exact S-matrix ([2]). 

The 1/N expansion of the SU{N) current form factor: We check the proposed 
exact form factor (l38l) in 1/A^ expansion. Fig. [5] shows the diagrams contributing to -F^M^f^ 
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Figure 5: Diagrams contributing to the form factor of the SU{N) current in the Gross- 
Neveu models up to order N'"^. 



in order and ^ which give 

Flk = €^\$f{q)ru{p)F\e) (45) 

(2vr)^ 



X <; v{q)Vi{k) 7± Vi{-k)u{p) - substr. }> + 0{N- 

+ 'jk — m 'jp — 'jk — m 



(46) 



The k integration can be performed using the propagators fHT]) and the vertices fHUl) . For 
convenience we write the total 4-point vertex function fj42|) which is a part of fH6l) as 



F = ^ A,,(fc)\/,(A;) ® V;(-A;) = F, + F, + Fy + F.^.^ 



with 



^ in 1 sinh( 

F. = — 1® 1 



iV cosh' (0/2) 
^ in ^ ^ 1 / sinh i 



iV sinh' (0/2) V 



ZTT 



rest 



%1X 1 

T7 7^ (7^^ ® 7^ (7^ - 2m) + 7^ (7/c + 2m) O 7^7^) 
iS k'^ 



where jk^jk = 'y^'jk^'j^'-fk+k'^'j'^^'j^ has been used. Correspondingly we decompose the 
form factor function F'^{6) in fH5l) (in order to avoid infra-red problems in the calculation) 
as 

F-\e) = 1 + (F^iO) + F^i9) + Fy{e) + FrUd)) + 0(iV-2). 
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The first contribution Ffj{9) is given by the 0(2A^)-Gross-Neveu form factor F^^{6) which 
has been calculated in 



F^e) = F^^(e) - 1 



1 / /-^ sinh^i^ /0cothi0-^cothi^ \A A 1 , ^ 



= ^ (^1 - 1^ - tanh i^) - 1^ f I + ^) - f I - ^) + ^ (i) 

^ 1 f § 
N ^sinh^ 

where 9 = in — 9 and -i/; (z) = (InF (z))'. Similarly we obtain 
1 / /""^ sinh^ i0 /0 coth i0 - ^ coth 



2iV \ Jo 02 + y cosh^ 10 - cosh^ 



and 



and therefore 



Frest{.&) = 



which agrees with the 1 /A^-expansion of the exact result for form factor of the the current 
derived in section [31 



5 Commutation rules 

In [31] commutation rules were derived for the Z{N) scaling Ising models. The results for 
the SU{N) Gross-Neveu model are similar, however, the proof is much more complicated 
because of the unusual crossing relations (120]) and (related to this) the Klein factors (ITOl) . 

Let I 9y^ with a = ((an, . . . , ain), • • • , (ctai, • • • , c^aro,)) be a state of a particles of 
rank ri, . . . , (1 < rj < N — 1) (or bound states of rj particles of rank 1). We define 
the charge of a state to be the sum of all ranks of the particles in the state 

a 

Qa = . 



6 CONCLUSIONS 



22 



The weight Wi{a) {1 < i < N) of the state a is equal to the number of ajk = i- Therefore 
the total charge of the state a is 

N 

Qa = J^Wj(a) 
i=l 

(see appendix [B]) . If a is a state for which the form factor F^{0 does not vanish we use 
and define the charge of the operator ip by 

N 

Qip = Qa mod N = wf mod 

with 0<Q^< N. 

Examples: For the energy momentum tensor T^'^(x) (which is a SU{N) scalar) the 
fundamental field ipaix) (which is a SU{N) vector) and the SU{N) current J^p^x) (which 
transforms as the adjoint representation) the weights and the charges are 

= (0,0,..., 0,0), Qt = 
= (1,0,...,0,0), Q^ = l 
w'^= (2,1,...,1,0), Qj = 0. 

Theorem 1 The equal time commutation rule of two fields (f){x) and ip{y) with charge 
and Q^, respectively, is (in general anyonic) 

4>{x)ij{y) = ^(y)0(x) exp {2me{x' - y^^ (1 - 1/N) Q^Q^) . (47) 

The proof of this theorem can be found in appendix [Bj 

6 Conclusions 

In this paper the general SU{N) form factor formula is constructed. As an application 
of this result exact SU {N) form factors for the field, the energy momentum tensor and 
the current operators are derived in detail. In the large N limit these form factors are 
compared with the 1/A^-expansion of the Gross- Neveu model and full agreement is found. 
The commutation rules of arbitrary fields are derived and in general anyonic behavior is 
found. We believe that our results may be relevant for the computation of correlation 
functions in fermionic ladders |3l]. In addition the series of the 1/A^-expansion of our 
exact form factors could hopefully help to understand the same series in QCD. 
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Appendix 

A Bound state form factors 

Proof of formulae (1291) — (1301) : Here we present a sketch of the proof for the bound 
state form factors formula. For simplicity several formulae will be written only up to 
constants, the normalization can be fixed at the end by the physical properties of the 
operator. The form factor formula for particles of rank 1 was proved in [13]. Applying 
the bound state procedure to this result we derive the formula for n = 1 bound state of 
rank N — 1, the general case n > 1 follows easily. 
The bound state intertwiner [28j is defined by 

t Res . ...z Res = if r(f) 

where the S-matrix S^{ip) exchanges all particles with rapidities (f = fN-i, ■ ■ ■ ,fi ^ f = 
ipi, . . . , ipN-i- It satisfies the bound state fusion equation 

r?^^£(^,^) = ^fi2(^,^)ri''^ (48) 

Lemma 2 The form factor for n particles a = ai, . . . , a„ of rank 1 and one bound state 
(p) = (p, . . . , pn-i) (with pi < ■ ■ ■ < pn-i) of rank N — \ may he written as 

Fa{p){0, uj) = yV2ij F^{9(p)T-^^ , with u = ^ _ ^ (v^i H h v^a^-i) 

for ipN~iN-2 = ■■■ = ^21 = irj, ^-e- ifj = u: + jir] - ire. 

Proof. We start with a form factor Fap,{6(p) for n + — 1 particles of rank 1 with 
rapidities 6 = 6i, . . . ,9n,^ = ^n-i, • • • , V^i and quantum numbers 
/Xi, . . . -iHn-i (for convenience we use for an inverse numbering). Applying iteratively 
the bound state fusion procedure (see e.g. [28|,|13]) we obtain the bound state form factor 

F^Ue,u){V2iY~''T^P^ =1 Res . . . ^ Res F^(^v^) (49) 

V / — ifiN-lN~2=iV </'21=«r? - — 

= F^{9ip)i Res ...i Res S^{^) 

— <^jV-lJV-2="? V>21=tV — ~ 

= i^^(^^)rf,)r(f) 
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where the form factor equation (i) Fg_^j^{6ip) = Fax{6ip)S^{ip) (see e.g. [13]) has been used. 
■ 

We start from the K-function K^^{9(p) for particles of rank 1 given by the general 
formula fl22|) where we replace ^ — > 6ip, {ip = ipi, . . . ,(pN~i) and integration variables 
z zy, {y = 1/1,.. . ,2/7V-2) 

KZ,ii.^) = j dz j dyh{9^,zy)p{eifi,zy) ^a{0^,zy). 

The state \E'a is a linear combination of the basic Bethe ansatz co-vectors f|TT]) (for n = 0) 

^aiOifi, zy) = L^^ {zy)^i{i^, zy) , with 1 < A 

where the L^^\zy) again satisfy a representation like the K'^y^{9(p). Iterating this nesting 
procedure we arrive at 

K'^xm = f dz^^ [ dy^'^ ... [ dz^"^''^ [ y^^'-'^hp^^^ 



where the functions h, and depend on the variables OLp,zy, {z = z^^\ . . . ,z^^ ^\ 
y = y}^\ . . .,y}^-'^\ = yf\ . . .,yN-i-i, (l = l,...,N - 2)). If we take the residues 
in fH9l) at = irj the pinching phenomenon (see [ISj and Fig. [1]) appears at y^^ = 

(pi, . . . ,^^^2 = This propagates to the higher level integrations such that we may 

replace y^^^ — > (^^''^ = Lfi, . . . , LfN^i^i which are related to cu hy ipj = u + jir] — in. The 
h-function (l23l) for the lowest level Bethe ansatz then takes the form (up to a constant) 

^, z^^\^^^) = 4>{i - z)4>{e - ^(^))0(^ - z)t{z)t{z - ^(1)) 
= /i(^, z)~^{l - -z)t{z- ^(1)) . 

Here and in the following we use the short notation 

n N-2 n 



m - ^'^) = n n - ^r)' ^u) = n 

i=l j = l l<*<i<"l 



T[Zi - Zj, 



et cetera, where the product is taken over all indices. The Bethe ansatz states defined by 
( ITT|) are related by 



where the bound state relation fl48p together with ^ and (jl]) has been used. These 
equations together imply 

u(a {^)\^f^^(a (i)^r^ + - ^) r ~fi{a) 
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The equations for bound state rapidity u — ipi — irj+in — (pN-i + irj — in and the relations 

b{z — <fj)(f>{z — Lpj) = — — z) and (f){ipi — z)/(t){z — </?jv-i) = —b{iTr — cu + z) have 
been used. Therefore we obtain the integral representation 



Kaip)iO,uj) = j dzh{9,z)p{9,uj,z)'^aip){0,uj,z) 



^a{p){0,i^,z) = L^pl^^{z,U})^^^~^{e,U},z) , with 1 < A, CTi = 1 < (72 < • • • < aN-1 



with the new K-, L- and p- functions given in terms of the old ones 



' (f)(^{2) + irj - z) ^^-^ ^ 
P{1:'^:^) = p{9(p, z(/^^ ) , {(pj = uj + jir] - in) . 
Correspondingly we obtain for a higher level Bethe ansatz / = 1, . . . , — 3 

and 

N-l 



i=N- 

such that 



and 



with / < 7i, ^1 = 1, . . . , Q = Z < Q+i < • • • < ^:Ar_i . 

Note that for / = iV — 3 the relation (p^cp^'-^'^^ +17] — ^^~^^^) = 1 holds because (p^^~^^ 
For the highest level I — N — 2 we have 
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such that 
and 



with ■ji = N, (Ti = 1, . . . , cr Ar_i = — 1 . 

Here 

L.-.M^x(.)^.(.O^r(i.^)r(i-l-|-). 

Finally we combine the minimal form factors in formula fl24l) for ipi^_iN_2 
(p2i = iv 



F^ipM,uj) = F^{e_^)Tl^ = F{e)G{e_-u)j — -K^ie^)^^^ 



(p{(p + iT] - 6) 

F{e)G{e-uj)K^(^pj{e,uj). 



The relations flTB]) for the minimal form factor function G for one particle of rank 1 and 
one of rank N — 1 and (fT7|) have been used. Therefore the final result is 

^Sp)(^'^) = / dz^'''...j dz^''-^'^h{e,u,^p'^{e,uj,^)^^^p){e_,u,£) 

N-2 riM-i 
Z=0 i=l 

P°{0,uj,^ = p'^ii^^gy) with = 

where p'^{Oip, zif) is the p-f unction for particles of rank 1 only. The complete Bethe ansatz 
state is 

I) = , u, ) . . . $«f ^ (.(^) , .(^) .(^) ) 

where is the highest weight bound state {rj) = {1,2, . . . , N — 1). 
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The energy momentum tensor: We apply the resuhs above to the example of the 
energy momentum tensor and prove (I34p . In this case n = n = 1, and the p-function is 
that of dSH) 

Lemma 3 The functions u) (for all I = 1, . . . , N — 3) are explicitly given as 

LpUi^^^) = ^M^^^K^ - ^) with 13 > I, (/i) = (1,2,...,/,*,...,*) 

L">(. - .) = e, r (- + —j r + —j (50) 

Proof. Again some equations are given up to unessential constants. We use induction, 
start with 

^ - .) ^ (-i)-r (1 + ^) r (i - 1 - ^ 

and then calculate iteratively for / = — 1, . . . , 2 the integrals 



'J Czuj 

'^iM^'^M^^''^ ^, m) = e^i^,)Sl\{z - u)Sf^^i{uj - u) 

= e^(^) [{N - I) 6^ph{z - u)d{uj -u) + 6fc{z - m)) 

where = 1 for {3 > I and else. Both integrals 

Ji = /" du(j){z - u)L^^'^ {to -u){N- l)h{z - u)d{uj - u) 



h = du(l){z — u)L^^\ijj — u)c{z — u) 
can be calculated by means of the formula 

^ dzT(a + —)T(b + —)T(c-—)T(d-—) 
V 27iiJ V 27iiJ V 2mJ \ 2niJ 



oo 



2 r (g + c) r (g + ci) r (6 + c) r (& + d) 

^ T{a + b + c + d) 



and yield the result (l50l) . ■ 
Finally we have to calculate 



K^^x^{9,u) = (e^^ + e"^) f dz ^{9 - z)L^'\uj ~ z)e'''esi,)^i%{9,u, z) 
esi^,)^%{9, ^, z) = e„(,) ((AT - 1) 6ib{9 - z)d{u - z) + 6>'c{9 - z)) 
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which yields the result using the formula 
for a + 6 + c + d = 0. 

B Commutation rules 

In this appendix we use the short notation for form factors i.e. matrix elements of the 
field ip{x) at a; = 

i^i%,0.) = F^i%,e^) = ^•°"*(^;3l^(o)l^)L". (51) 

To proof the general commutation rules of fields (1471) we have to consider SU{N) sum 
rules and general crossing relations. 

B.l SU{N) sum rules 

Particles and anti-particles: Let (a) = (ai, . . . , a^), (1 < ai < ■ ■ ■ < < N) a 
particle of rank (and charge) r. The corresponding anti-particle is (a) = (ai, . . . , d]\f_r), 
(1 < ai < ■ ■ ■ < aN-r < N) (of rank N — r) such that the union of the set of indices 
satisfies {ai, . . . , a^} U {di, . . . , dN^r} = {!,•••, N}. Therefore 

r N-r N ^ 

Y^a, + J2dk = J2k = -N{N + l). (52) 

k=l k=l k=l 

Charges: Let a = ((cm, . . . , air J, • • • , (ooi, . . . , ctara)) be a state of a particles of 
rank ri, . . . , (1 < r j < — 1) (or bound states of rj particles of rank 1). We define 
the charge of a state as the sum of all ranks of the particles in the state a 

a 

Qa = T'j . 

The charge of anti-particles (bound states) we define as 

a 

Q^ = 5^r,(iV-l) = (iV-l)Q„ 

Qa ~l~ Qa -^Qa • 
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Weights: Let a — ((ccn, . . . , ccin), . . . , (cKqi, . . . , Oiara)) be a state of a particles. The 
weight Wi{a) of the state a is equal to the number of ajk — i, {1 < i < N) 

j=i k=i 

Therefore the total charge of the state a is 

a a i'j a ''j N N 

j=l j=l k=\ j=l k=l i=l i=l 

Similarly, we consider 7 = ((711, . . . , 71^^, (7-^1, ... , 7^^^)). 
Sum rules: Because of SU (N) invariance 

(or ■ipa'y 7^ 0) implies for the weights 

w{a) = w(7) +w^ + L{l,...,l), L e Z 

where w'^ is the weight vector of the operator t/j and (1, . . . , 1) are weights of a state in 
the vacuum sector. Therefore 

N N 

i=l i=l 
The charge of the operator t/j is defined by 

= {Qa - Q-y) mod iV, Q<Q^<N (53) 

AT 

= wf mod N . 
1=1 

For a particle (ttji, . . . , ocjrj) of rank Vj we use the short notation (aj) = (oji, . . . , ctjrj) 
and ctj = X^fcLi Oijk and correspondingly, jj — Yll^=i Ijk- Then SU{N) invariance implies 

a 7 

With R^^J2i^t-liN + l)lj2'^f-QA (54) 

i=l \j=l / 

which can be straightforwardly proved using the above definitions. 
Examples: 

rpnu . ^ (o,0,...,0) Qt^O Rt^O 

t/ja-. = (1, 0, . . . , 0) = 1 i?^ = 1 
j^u . ^ (2,1,..., 1,0) Qj = Rj^l-N 
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B.2 Crossing 

B.2.1 A partial S-matrix 

Definition 4 Let 9_p = (^^7r(i); • • • ? ^7r(a)) be a permutation of 6^ = {6i, . . . , 6^)- Then 

Sa{0_i3;9_^) is the matrix representation of the permutation group Sa generated by the 
simple transpositions aij : i ^ j for any pair of nearest neighbor indices I < i,j = i + l < 
a ai§ 

c^ij S{6ij) 

Because of the Yang-Baxter relation and unitarity of the S-matrix the representation is 
well defined. We will also use the notation 

if IT is that permutation which reorders the array 9^ such that it coincides with the combined 
arrays of and 9_y^. 

As an example consider the case = (6'i, 6*2, 6*3, 6*4), 9^ = (6*2, 6*3) and 9_x = (6*1, 6*4) 

s^-{e,eM; 6^626,6,) = s^^'^i{en)Si^i^i{eu)s'j^ 

/il ^2 Ai Ai 

102 1^3 1^1 1^4 ^1 ^2 Ai A2 

I g I - V< I 

1^1 1^2 1^3 ^4 ttl 0(2 Ct3 «4 
«! a2 tts a4 

If the permutation inverts the rapidities completely S§; [6_ij; 6^) = S§_{6^) is the full S- 
matrix. 

B.2. 2 Crossing for SU(N) 

As was argued by Swieca et al. [8| the particles of the chiral SU (A) Gross-Neveu model 
posses anyonic statistics and due to the unusual crossing property of the S-matrix, Klein 
factors are needed. The crossing relations of the form factors for normal fields and particles 
were derived in by means of LSZ-assumptions and maximal analyticity. They have 
to be modified for the chiral SU{N) Gross-Neveu model. 
We propose crossing relations 

i^M^;L) = 4) E (55) 
= E C,.)^i^(^7;^4<^)C^%i^#(M,;^) (56) 



^Note that this definition is quite analogous to that of representations of the braid group by means of 
spectral parameter independent R-matrices. 
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which, compared to the formulae in [2H], are modified by the factors (f ^ s and ^ 



-y pi^{N-l)^Q^{Qr,+N)^{R^Q^-Q^Y:ij) 
^(7,a,r,) - P{^,af ^ 

^ „iT(Af-l)(|Qr,(Q„+Af)+Q.Q^)^(fi^Q<.-Qv'E7j) 

'^(7,a,r)) 

= ('_1^(^-l+(l-W)(Q«+'^7-<^^))Q7 

" (7, a) ' 
^'A ^ i7r(l-l/Af)Q^Q^ 
{7) 

with 7j = ^N{N + 1) — 7j, due to (152]) . The sign factor p^^^-^ and the statistics factor 
(T^j were introduced in [13j. The charge of the operator ijj and the number are 
defined in fl53D and fl54D. 



B.3 Commutation rules 

In [21] commutation rules were derived for the Z{N) scaling Ising models. The results for 
the SU{N) Gross-Neveu model are very similar, however the proof is more complicated 
because of the unusual crossing relations and the presence of the Klein factors. 

Theorem 5 The equal time commutation rule of two fields (f){x) and ip{y) with charge 
and Q^, respectively, is (in general anyonic) 

(j){x)tlj{y) = ^{y)(j){x) exp {27ne{x' - y')l (1 - 1/N) Q^Q^) . 

Proof. We consider an arbitrary matrix element of products of fields 

immfM^e^) = ^'"''\o:,\<f>{xmy)\e^y^. 

Inserting a complete set of intermediate states | )™ we obtain 

JP'x-iP^y]_ / a \j.-l(n n ^ -iP-,{x-y) 



where Pa = the total momentum of the state | )^ etc. and = 111=1 / i^- Einstein 

summation convention over all sets 7 is assumed. We also define 7! = H^Li ^r^- where 
is the number of particles of rank r in 7. We apply the general crossing formulae (15511561) . 
Strictly speaking, we apply the second version (156|) of the crossing formula to the matrix 
element of ch 



where p = [pp, . . . , pi) with p = antiparticle of p and 6^^ — iti- means that all rapidities 
taken the values 6'' — z (tt — e). The matrix is defined by flSTj) with O = 1 the 
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unit operator. Summation is over all decompositions of the sets of rapidities O^/j and 9^. 
To the matrix element of ip we apply the first version of the crossing formula fl55l) 



We insert fl55l) and fl56|l in fl57|) and use the product formula S^~{6^6^]9^)Suri{0^]0.y0ri) = 
S%{9^0^] 9_^9_j^) . Let us first assume that the sets rapidities in the initial state 6^ and 
the ones of the final state 9^^ have no common elements. This also implies H ^ = 0. 
Then we may use (ii) to get {9^9^; 9 ,^9^) = 5^ and then we can perform the 9,^- and 
^-integrations. The remaining ^-integration variables are 9^ = 9_^ 9_^- Then we may 
write for the sets of particles £ = (ilj, t] = urr and 7 = ^urr and similar for rapidities and 
momenta. Equation (1571) simplifies to 

X 5i-(^)e<^^^'')^-^(^-^^)^ (58) 

where 

X C^C^C^7/>^(^; + ^7r„,|^ + z7r_,^Je-^-^-("-^) . (59) 
Similarly, if we apply for the operator product ip{ii)(t){x) and 

use the second crossing formula to the matrix element of (f) 
and the first one to the matrix element of ip 

Similarly we obtain (with x = V = ^ = iiijJt) equation fl58|l where replaced 
by 

X C^C^C^V'f^(^; + z7r_,£^,^,)e^^-(^-^) (60) 
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which means that^only aJ'^^^^Cl,a,v)Hfi,'y,P) replaced by (Tlp^s)Cw,5,rM5,a,^) and the inte- 
gration variables 0^ hj 9^ — 27r_, i.e. P^j by —Puj. 

If there were no bound states, there would be no singularities in the physical strip 
and we could shift in the matrix element of %l){y)(f){x) fl58l) with fl60|l for equal times and 

< the integration variables by 9i 9i + i7i_. Note that the factor e*^'^*^^"^) decreases 
for < Re6'j < vr if < y^. Because P^ —P^ (if ^ £i, ~ i^^-) we get the matrix 
element of (p{x)ip{y) fl58|) with fl59|) up to the factor 

(7.«)^(7,«.^)^(/3,7,P) ^ g-27rii(l-i)Q^Q^ 
"(/3,5)''(/3,5,r)'S((5,a,<^) 

This equality follows after a long and cumbersome but straightforward calculation. In 
pi] was shown that we obtain the same result if there are bound states. ■ 
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